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Abstract 
> _ 

r^ I Shell structure in the single particle spectrum of deformed har- 

monic oscillator potentials when a term proportional to (L)^ is added 
CSJ ■ is analyzed for a large particle number. In the case study presented 

J~^ [ here it is argued that, in view of the chaotic nature of the problem, 

Qv^ I a thorough understanding of the classical situation provides essential 

guidance in tackling the corresponding quantum mechanical problem. 
A scaling law which gives a dividing line between regular and chaotic 
behavior in terms of energy, deformation and strength of the (L)^ 
'^ ■ term has been found. According to this law, shell structure survives 

fi , for higher particle numbers only with lesser deformation. 

O 

> 

'O : PACS numbers: 05.45. +b, 36.40.Qv, 71.24.+q, 21.10.-k 

1 Introduction 

Finite Fermi systems give rise to shell structure as is well known from the 
existence of magic numbers in atoms and nuclei. Nowadays, shell effects are 
also established in metallic clusters [|l], |^, ^. The first experimental evidence 
in alkali metallic clusters 0] , and independently the theoretical prediction of 
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magic numbers within the framework of the spherical jelhum model p, ^], 
have been major hallmarks in this field. 

The valence electrons in metallic clusters are strongly delocalized. They 
play an essential role in the understanding of many observables and they can 
be treated by the mean field approach. The jellium model works very effec- 
tively for monovalent metallic clusters like sodium and potassium. This mean 
field model describes a single valence electron that interacts with the average 
potential which is generated by the other electrons and ions. The structure 
of the ions is neglected, they are assumed to form a constant positive back- 
ground or jellium density within a finite volume. While the spherical jellium 
model has been able to explain general features of enhanced abundances of 
valence electrons in spherical clusters 0-0, it has been recognized that the 
finer structure in the mass spectra between magic numbers can be explained 
via symmetry breaking mechanisms [0]-||T2| similar to the situation in nu- 
clear physics |T3|. A spherically symmetric mean field leads to very strong 



shell effects which is manifested in the stability of the noble gases, magic 
nuclei and magic metallic clusters. When a spherical shell is only partially 
filled, a spontaneous breaking of spherical symmetry can give rise to a de- 
formed equilibrium shape, if it is associated with an energy gain. For clusters 
with incompletely filled shell, it has been suggested, in analogy with atomic 
nuclei, to consider a quadrupole deformation which leads to the lifting of the 
degeneracies in the spherical model. In fact, this idea was implemented by 



Clemenger [|1^ who applied the modified harmonic oscillator model (MHO) 
to metallic clusters where the model played an important role in the under- 
standing of various experimental data (see for review |]1], 0). The MHO with 



a spin-orbit term was introduced by Nilsson [Q in nuclear physics, and it 
is a very effective tool in the analysis of ground state deformations and sin- 



gle particle excitations for many nuclei |T3|, |T5|. The same model without 
the spin-orbit term seems to work successfully in reproducing the shapes of 
small clusters which have been calculated by the more sophisticated meth- 
ods of quantum chemistry ^. Based on these results it is therefore accepted 
that cluster deformations exist. This is in fact confirmed at least for clusters 
with A < 40 either by the Clemenger-Nilsson (CN) model or by a self- 
consistent Kohn-Sham density-functional method |16| (KS) with deformed 
jellium backgrounds ||l3-[|20l . 



Nowadays there is experimental evidence of shells for metallic clusters 
with more than 1000 atoms [^-0,|21|, ^. The shell ground state proper- 



ties of spherical metallic clusters containing a few thousand valence electrons 
can be described within the KS approach using the spherical jellium model 



[ p^ . Solving the KS equations is quite difficult when deformation effects are 
taken into account, even for modern supercomputers. It turns out that phe- 
nomenological potentials used traditionally in nuclear physics serve a purpose 
similar to those obtained within the KS approach if the relevant parameters 
are adjusted appropriately. Typical potentials are the Woods-Saxon (and 
its various modifications) (WS) and the MHO potential, i.e. the CN model. 
The main idea of the approach exploiting the WS and the MHO is to use the 
single particle spectrum generated by the potential to calculate shell energies 
only, because neither WS nor CN models allow for a description of binding 



energies for neutral or charged clusters. Following Strutinsky's idea [^, the 



energy of an interacting Fermi system is divided into a smooth part, denoted 
by E, and an oscillating part which is the shell correction energy denoted by 
6E, i.e. E = E + 6E. The quantity E varies slowly with particle number and 
the deformation of the system. In nuclear physics and for metallic clusters it 
is usually replaced by the liquid drop energy which carries the information 



about the bulk properties of the system. The shell correction energy 5E con- 
tains all the oscillations originating from the bunching of the energy levels. 
It is obtained by summing the single particle energies of a phenomenological 
shell model potential and subtracting the average smooth part of the total 
energy 0, |l|, ||]. 

Using the spherical WS potential for clusters with thousand valence elec- 
trons Nishioka et al ^5| predicted the super-shell phenomenon where the 



regular oscillations in the level density which reflect the main shells are mod- 
ulated by an amplitude which oscillates with a lower frequency. Such analysis 
led to the experimental discovery of the most stable large spherical metallic 
clusters ||2l|, ^ |^. Employing the Strutinsky shell correction method, a sys- 
tematic study of cluster equilibrium shapes has been performed for the WS 
potential |]1T], ^ up to 300 atoms and for the MHO potential for clusters 
with particle numbers up to 300 in [^, ^ and up to 800 in |T^. Effects of 



deformation are clearly seen in these calculations. For clusters with particle 
number of up to a few thousand it is natural to use these potentials for the 
analysis of the existence of deformation. The problem is still being discussed 
in the ongoing literature (see, for example, discussions by S.Bj0rnholm in [^ 
and H). 

The equilibrium deformation is ultimately related to the behavior of the 
shell correction energy 5E of the quantum Hamiltonian. According to the 



periodic orbit theory |Q-|^ (see also |jT3|, |3^) the frequencies in the oscilla- 
tions of quantum single particle spectra are determined by the corresponding 
periods of classical closed periodic orbits. The gross shell structure is due 
to the periodic orbits with a shortest period, whereas the finer details are 
determined by contributions of longer orbits. This result has been confirmed 
for spherical potentials used for large clusters PH pil BH BB]. It suggests that 



the semi-classical approach can provide deeper insight also into the problem 
of large deformed clusters by exploiting the tools which have been developed 
in the study of classical integrable and non-integrable systems. 

One of the first attempts to analyze the interplay between shell effects 
and classical chaotic motion has been done in ||3^ for an ellipsoidally de- 



formed version of a diffuse nuclear average potential called the Buck-Pilt 
potential. If the classical problem is non-integrable and displays chaotic be- 
havior, the shell structure of the corresponding quantum spectrum is affected 
depending on the degree of chaos. More recently these effects were demon- 
strated for the quadrupole + octupole |^, quadrupole + hexadecapole |^ 



and for quadrupole + octupole + hexadecapole ||4^ axial deformed harmonic 
oscillator potentials. 

Initiated by these latter investigations we make use of the connection be- 
tween ordered motion in the classical problem and shell structure phenomena 
in the corresponding quantum spectra for phenomenological potentials used 
for clusters. In practical applications the study of the more realistic de- 
formed WS potential is complicated due to the effective occurrence of many 
multipoles; this is the subject of a future paper. Since the MHO is being 
used widely for the interpretation of data |jl|, |12|, ^ ^, Q we expand in 



the present paper the preliminary classical analysis of the MHO potential 



4ip and demonstrate its connection with quantum-mechanical single parti- 
cle orbitals. One particular result of the classical analysis is a kind of scaling 
relationship between energy range and the strength A of the (L)^ term in 
the MHO potential. For the corresponding quantum mechanical problem 
this translates into the finding that pronounced shell structure can occur for 
increasing particle numbers only for an accordingly decreasing strength A 
and/or a decreasing quadrupole deformation. 



2 General properties of the model 

The level ordering in the WS potential falls between the soft-surface harmonic 
oscillator (HO) and the hard-surface square well. A similar level ordering is 
obtained by modifying the HO. The MHO reads in cylindrical coordinates 

H = ^iPl+Pl + f ) + ^(^^ + ^) - ^^-iLf (1) 

where L is the angular momentum in units of ft and A > 0. As the model 
was originally designed within a pure quantum mechanical context [jl3|, |1^] we 
conveniently use the energy unit huo for the (L)^-term. We consider only axial 
quadrupole deformations. For 6 = 1 all three components of L are conserved 
because of spherical symmetry; the model becomes then integrable. For 6 > 1 
and 6 < 1, which is referred to as prolate and oblate quadrupole deformation, 
respectively, the problem is non-integrable. The axial symmetry makes the 
problem a two degrees of freedom system with L^, denoted classically by p^, 
being a constant of motion. 

The classical phase space is unbounded for given energy. The subtrac- 
tion of a term proportional to the square of the angular momentum allows 
in principle arbitrarily large values of positions and momenta. This is easily 
recognized for 6 = 1 where an arbitrarily large value of the then conserved 
quantity (L)^ can be given as one of the initial condition; using correspond- 
ingly large values of positions and/or momenta, an arbitrary value of the 
energy can be attained. The corresponding quantum mechanical levels are 
given by the expression En ~ hujl{l — XI) for rir = {N = 2nr + I) for 
sufficiently large /. For / > 1/A the energy levels are therefore negative. This 
unboundedness has been dealt with tentatively by adding to the Hamiltonian 
a term proportional to (A^|L^|A^). The expectation value is taken with shell 
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model states and can be expressed in terms of shell numbers |T3|, |T5[, if such 
quantum numbers are well defined which applies to the regular region only. 
For & 7^ 1, when (L)^ is no longer conserved, its time variations extent over a 
fairly large range of values for generic orbits; depending on initial conditions 
such values have no bound in principle. We enter the chaotic region and 
the corresponding quantum mechanical states become mixtures of angular 
momentum states. As long as the range of such admixtures obeys the in- 
equality / < /crit with A • /crit = 1 foi a givcu value of A, the states are not 
affected by the unboundedness of the Hamiltonian (|I]). It is discussed below 
that, for 6 > 1 and the range of A- values considered, such states give rise to 
the chaotic nature of the spectrum. In contrast, states with large values of 
/ exceeding /crit are basically regular and do not mix with the former group 
of states. The latter states give rise to negative energies but are irrelevant 
for our discussion as they are essentially orthogonal to the former relevant 
group. 

Since the model is non-integrable for 6 7^ 1 and A 7^ there is no obvious 
basis for the corresponding quantum mechanical problem. If A <^ 1 and 
b > 1 the adequate basis makes use of the asymptotic quantum numbers 
Q which diagonalize the Hamiltonian for A = (see below). While this 



approach facilitates the technical problem, even this representation requires 
great numerical effort to tackle the problem. In accordance with the classical 
problem, the quantum spectrum shows all characteristics of chaotic behavior 
for large systems. 

We have chosen the upper limit 6 = 2, which is denoted super-deformation 
for nuclei [|1^. Super-deformation seems to exist in metallic clusters, at 
least for small particle numbers (see, for example [1|, ^ |19|, |2y, ^, ^ p8|). 
We focus our attention only on 1 < 6 < 2, since b < 1 (oblate case) can 
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be obtained from the results of the former case by interchanging q and z 
and suitably replacing uj hy huj. This is in stark contrast with the situation 
considered in |^, |5^, ^. There the oblate case is chaotic and shell structure 
is essentially destroyed once higher multipoles are added while shell structure 



can prevail in the prolate case. Note, that the potential used in ||3^ contains 
effectively higher multipoles. 

3 Classical treatment 

For most of the results presented here it suffices to consider p^ = in Eq.(|l]). 
Larger values of p<^ are discussed at the end of this section. 

We first consider qualitatively the case of spherical symmetry, i.e. 6=1. 
The problem reduces to one degree of freedom since now pg = mOr'^ = 
zpg — gpz is conserved. Here r and 6 refer to the usual polar coordinates. 
Closed orbits occur if the radial and angular frequencies are commensurate 
(see, for example |0). Rewriting Eq.(|I]) as 

H = {pI+ pl/r'^)/{2m) + muj'^r'^/2 - Xupl/h (2) 

we obtain closed orbits if 

t/s = 1/2 - Xpe/h (3) 

with integers t, s. For instance, when t/s = 1/3, 1/4, . . . the trajectory forms 
essentially a triangle, a square etc. in the g — 2— plane; for t/s = 2/5 we get 
the five star, and so forth. The precise shapes of these geometrical figures 
depend on the magnitude of A in that small values of A produce polygons with 
rounded corners while larger values yield loops at the corners as is illustrated 
in Fig.(l). The appropriate scaling of A is given by the kinematical constraint 



between the energy and the angular momentum which reads 

E + Xupyh>uj\pe\ (4) 

Here we stress again that for positive values of A, as used in actual applica- 
tions, the value of the angular momentum pe is not limited in its absolute 
value for given energy; in particular, if 4A|-E| > Tiuj, arbitrarily large positive 
or negative values may be assumed by po. 

For the deformed system, the orbits are Lissajous figures when A = 0. 
Closed orbits are obtained for rational values of h. Our interest is directed 
towards short periodic orbits. Fig. (2a) displays the phase space structure 
at f) = where the surfaces of section are taken. Note that the whole 
region to the left and right, i.e. outside the lines z = ±Jfi/{2Xmuj), is 
accessible as long as the two lines do not intersect with the ellipse which 
forms the other part of accessible phase space. If A is sufficiently large to allow 
intersection of the lines with the ellipse, the phase space becomes connected. 
In this case the part of the ellipse which is to the left and right of the lines 
is inaccessible. In Fig. (2b) surfaces of sections for A = 0.03 are given for 
four different orbits. With the choice of energy E = 10 and the frequency 
u = 1.5 the two lines are outside the ellipse. (For convenience we use units 
for which m = 1 and hu = 1). The same pattern is obtained if E and A 
are rescaled such that EX =const (see discussion below). The diagram of 
Fig. (2b) refers to 6 = 2, but the pattern as described prevails for 1 < 6 < 2. 
Note that, in accordance with the estimates below, the onset of chaos is 
discernible just at the periphery of the ellipse, which is the region where Eg 
(the energy associated with the motion in the ^-variable) is small. Typical 
trajectories of periodic orbits are displayed in Figs. (3); Fig. (3a) corresponds 
to the black dots on the far right and left of the phase diagram. Fig. (3b) is 
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the orbit sitting in the centre of the four islands, and Figs. (3c) and (3d) are 
represented by crosses in Fig. (2b). It was observed in |^ that there is a 
remarkable similarity between short periodic orbits discussed above and the 
ones occurring in a WS potential with increasing quadrupole deformation. 
In fact, the short periodic orbits are of the same geometrical shape and 
occur in similar regions of the surfaces of section. This similarity could 
indicate that shell structure effects may be similar for the two potentials. 
A more careful analysis taking into account stability and degeneracies of 
corresponding orbits [jl3|, ^ goes beyond the scope of this paper. 

For larger values of A hard chaos takes over quickly within the ellipse, 
in particular when the two lines enter the ellipse. We make two important 
observations: (1) regular motion prevails far outside the ellipse which is the 
whole area to the left of the left line and to the right of the right line. 
These orbits may attain very large values of q and z; also the variation of 
the angular momentum pg is unlimited in principle. (2) The variation of pe 
ranges typically between —50 and +50 for generic orbits inside the ellipse for 
our choice of parameters. In Figs. (4) we illustrate typical surfaces of section 
and trajectories for the two cases. The orbit on the bottom left of Fig. (4) is 
representative for a regular orbit and its surface of section produces the kind 
of tori lying completely outside the ellipse as illustrated on the top left of 
Fig. (4); the irregular trajectory on the bottom right of Fig. (4) produces the 
irregular dots inside and outside the ellipse as shown on the top right; note 
the proximity to the periphery of the ellipse of the dots outside the ellipse. 

The characteristics of these two types of orbits are significant when con- 
sidering a corresponding quantum mechanical calculation. The nearly regular 
motion for very large values of pq is expected to lead to little mixing of the 
class of states which is characterized by very large (/ > 60) and fairly sharp 
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values of the angular momentum. Also the quantum shell number is still well 
defined for these states. These are the states which make the quantum me- 
chanical problem unbounded from below. Owing to their negligible mixing 
to the states discussed in the following they are virtually orthogonal upon 
the states in the chaotic region. Their classically regular nature is refiected 
by the de facto existence of good quantum numbers. The occurrence of chaos 
is therefore not the result of unbounded motion, neither classically nor quan- 
tum mechanically. However, the other class of states is expected to have 
not only substantial expectation values of the operator (L)^ (corresponding 
to / ~ 30) but in particular a variance of similar magnitude. Consequently, 
shell structure is expected to be quite weak and the concept of shell quantum 
number becomes ill defined. There is of course a gradual transition between 
the two classes. The values given for / refer to A < 0.02. We note that these 
physically relevant states obey the criterium / ■ A < 1. The results of the 
following section fully confirm such expectations. 

We conclude the discussion of the classical motion by the results of a 
perturbative treatment ESI which was successfully employed in HDI. When 



6 > 1 we average over the angle 6g after rewriting the Hamiltonian in terms 
of action and angle variables of the unperturbed problem. This is sensible as 
long as the variation of 6g of the exact motion is faster than that of 6^, which 
is ensured if E^^ E^, where Eg and E^ are the respective energies residing 
in the q— and z— motion. In principle, these energies are not constants 
of motion, but the two degrees of freedom become uncoupled within the 
approximation. In fact, we obtain after averaging of Eq.(|l]) the effective 
Hamiltonian H^^ = Hg + Hz with 

— ^ 

2m 2 



Pg rnw\2 
-"e = 7717 ^ 1^& 
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(5) 



^2m 2 ^ 



where 

a = 1-2A^ 



(6) 



^' = 4h-''bi" 



To ensure a real effective frequency cu^ in Eq.(|^) we read off an inequality 
which must be obeyed to render the approximation valid. In view oi E ^ Eg 
it may be written as 

AS < W- P' 

We note that the exact solutions are perfectly well represented by the motion 
governed by H^s if the inequalities are observed. We return to the important 
inequahty (|^) in the next section. 

So far we have concentrated our discussion upon p^^ = which suffices 
for low values oi p^. For larger values the pattern changes, however, in that 
chaotic motion is suppressed to an increasing extent for increasing values of 
p^. Stable periodic orbits occur for large values of p^ at A-values for which 
the situation is fully chaotic at the lower p^-range. An example is illustrated 
in Fig. (5). These orbits will of course prevail in the corresponding quantum 
mechanical spectrum, when the full spectrum is considered. We discuss in 
the next section that such levels can lead to an overestimation of shell effects 
for the total energy of large clusters. 
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4 Quantum mechanical treatment 

As it was mentioned above for deformed system the obvious basis uses the 



asymptotic quantum numbers [jl4|, ^ which diagonahze the Hamihonian for 



A = 0. The boson operators are a^ and aj for the z-motion while for the two 
dimensional motion perpendicular to the 2;-direction we use 

A± = —=[a^^%ay) (8) 

A = ^(4±mt). (9) 

The conserved 2;-component of the angular momentum reads in this basis 

Lz = n^ — n_ (10) 

where rij- are the occupation numbers referring to the operators {A±yA±. 
For a fixed eigenvalue of L^, denoted by A, we effectively have two quantum 
numbers defining the basis, i.e. n^ and ra+. This reflects the two degrees of 
freedom of the corresponding classical situation. 
The unperturbed (A = 0) spectrum is given by 

2t7 -I- 1 
K.,n+ = nu{2n+ + 1 - A + -^^) (11) 

and the square of the angular momentum operator reads 

(Lr = l{l + b){{2nz + l){2n^ + l-A)-{al + {air){D + D^)) 

+ lil - h){{al + (4)2)(2n+ + 1 - A) - (Z^ + D^){2nz + 1)) 
+ {al-{a\f){D-D^)-l + X' (12) 

where we used the combination D = A+A^. 
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It is now straightforward to set up the matrix of the full Hamiltonian. 
Parity conservation is reflected by the fact that the problem reduces into 
even or odd occupation numbers n^. The occupation numbers n_|_ start with 
the positive value of A. We stress that the general form for (L)^ as given 
by Eq.(|l^) mixes states with AA^ = 0,2,4 for 6 7^ 1. As a consequence, 
the calculation cannot be restricted to a particular shell number A^ for the 
deformed case as it is usually done in nuclear physics ||13|, |l^. In particular. 



the wave functions will be superpositions over a range of shell numbers where 
the degree of the mixing increases with increasing deformation and/or energy. 
From Eq.(p!ID we conclude that, if one decides to truncate at n^ = N^, 
the range of n^ must extent up to A^^ = [2b ■ N^] where [...] denotes the 
integer part of a given number. With this choice it is ensured that the unper- 
turbed basis is consistently represented up to -Econs = (2A^+ + 1 — A + 1/(26)) 
(in units of hu = 1), in other words, all possible degeneracies are taken 
into account up to Econs', for 6 = 1 these degeneracies are just the pos- 
sible angular momentum values. Note that the range between -Econs smd 
Ejaax = (2A^+ + 1 — A + {2Nz + l)/(26)) cannot be omitted for reasons of 
consistency, since n_|_ and n^ must run independently to their respective max- 
imal values. These considerations are crucial in the present context as the 
mixing of low lying levels with high angular momenta is essential. It is in- 
structive to discuss first the mixing of angular momenta for 6 = 2 and A = 
where the eigenstates are exactly available. In Table (1) we give the aver- 



age L^ = {n\{Ly\n) and the variance AL^ = y {n\{L)'^\n) — {L'^y for some 
selected positive parity levels for A = (recall A = and huj = 1). Note in 
particular the large values of the variances which often exceed those of the 
averages. The switching on of the (L)^ term not only lifts the degeneracies 
but brings about a strong mixing among all levels. This becomes obvious in 
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Table (2) where we display corresponding quantities for A = 0.01. The ap- 
preciable mixing of high angular momenta signals interaction of many levels 
on a large scale which is also confirmed by the eigenvectors. For the numbers 
chosen they pick up typically twenty and more components larger than 0.1 
from rather distant unperturbed states, the eigenstates become strongly de- 
localized. This mechanism requires a considerable extension of the basis in 
mean field calculations even for relatively small strongly deformed clusters. 
The strong mixing between levels, irrespective of the basis chosen, is re- 
lated to the onset of chaos ^^ , which is confirmed by the nearest neighbor 



analysis of the spectrum. A particular sample is illustrated in Fig. (6); a sam- 
ple average would yield an almost perfect match. This strong result cannot 
be affected by a volume conservation condition lo^lOz = const., which is used 
for the calculation of Potential Energy Surfaces in order to find global and 
local minima, which yield equilibrium shapes and shape isomers, respectively 
(see e.g. [^, |2^). In fact, while the effect of this constraint may change the 
level ordering once the (L)^ term is included, it is observed (e.g. [|T^) that the 
ellipsoidal shape at large distortion becomes unfavorable for any combination 
of filled single particle orbitals. In addition, the increased level density at 
higher energy enhances the interaction of the single particle states. There- 
fore, there is little scope to find pronounced shell effects for a system with 
particle numbers that fill shells beyond A^ = 4 in the unperturbed basis and 
b > 2. Likewise, for A^ > 20, even a small deformation such as 6 ~ 1.1 does 
not support shell structure. 

The levels listed in Table (2) are selected according to the strength of 
their eigenvectors being essentially concentrated in the upper fifth of their 
total length. We have chosen the truncation of n^ at A^+ = 22 which gives a 
1035 • 1035 total matrix size. The states listed have half of their total norm 
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exhausted within the first 207 components. The spectrum of the states so 
selected is displayed in Fig. (7a) as a function of A. In the top right corner 
we discern a distinctly irregular pattern which is caused by the arbitrary 
criterion applied to select the states. In fact, in this region the eigenvectors 
depend very sensitively on the strength A, since this is the region where 
level repulsions occur on a large scale. Therefore, states from high lying 
unperturbed levels will, owing to their coupling to lower lying states, pick up 
considerable strength (0.5 for our choice) in their first 207 components for a 
certain window of A-values while other states will undergo a similar change 
for a different window of A-values. In Fig. (7c) we display a section of the 
same states without connecting the data points. Here it is clearly seen that 
a few states maintain their identity for increasing A while the majority of 
states is submerged in the morass of a chaotic spectrum. To make this point 
even clearer, a stricter criterion is used in Fig.(7d) where data are given for 
state vectors which exhaust 80% of their norm in the first 207 components. 
Obviously, only a few levels survive, in fact about just the ones which mix 
weakly. Notice the fairly sharp dividing curve between regular and irregular 
behavior in Fig. (7a) which, according to the previous section, should be given 
by the scaling law EX =const. 

When no selection is made there are many more levels with comparable 
and even lower energies for A > 0.005 as is illustrated in Fig. (7b) where 
the complete spectrum is displayed. These additional states mix to a lesser 
extent with the ones discussed above. The additional levels come from higher 
lying unperturbed levels as seen in Fig. (7b) and have even higher angular 
momentum components. These states are nearly orthogonal to the states 
belonging to the chaotic spectrum; they correspond to the states discussed 
in the previous section and form the space giving rise to the negative part 
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of the energy spectrum. As a function of A the energies come down owing 
to the negative (L)^-term; accordingly, the corresponding eigenvectors have 
their bulk components at higher labels. They are the quantum mechanical 
counterpart of the type of classical orbits displayed on the left column of 
Fig. (4). Traditionally these states are pushed up in their energy by the ad 
hoc adding of the mean value of L^ which is, for small deformation, a simple 



expression in the shell number |T^, ^. This remedy does not affect the 



eigenvectors, it would change the spectrum but not its statistical properties. 
The basic reason for this to hold is the fact that the states whose energies 
tend to minus infinity do not couple with the physically relevant states as the 
irrelevant states have very large angular momentum while the relevant states 
have a mixture of only moderate angular momenta. Since the concept of 
shell numbers breaks down, the mean value of L^ can no longer be expressed 
in terms of a shell number and adding it to the spectrum becomes doubtful. 
On the basis of our findings we rather ignore the irrelevant states. 

For completeness we report corresponding results for 6 = 1.1. In Table 
(3) we present results for A = 0.1. To demonstrate the only slightly bro- 
ken symmetry for the near spherical case this value is chosen much larger 
than the one usually applied in either nuclear physics or metallic clusters. 
As a result, irrelevant states with large mean values of angular momentum 
(/ > /crit = 1/A) have negative energies. All states are fairly good eigenstates 
of angular momentum with much smaller spreading AL? as is expected for 
this mild deviation from unity of h. In most cases the expectation value of 
(L)^ is very close to an eigenvalue l{l + 1), we have indicated the correspond- 
ing /-values in brackets. The eigenvectors are fairly well localized, there is 
only limited mixing. A near Poisson distribution is found for the nearest 
neighbor distribution indicating the near-integrability of the 6 = 1.1 case. 
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Note, however, that these statements are hmited to about the first three 
hundred levels. 

Based upon the inequality E^ ^ Eg we obtain from Eq.(|^ the border 
line EX ~ 1/(26^) between regular motion and expected chaos. A similar 
dividing line should be expected for the quantum mechanical problem. It 
appears that the relation EX ~ 1/(46^) can serve as an empirical dividing 
line between order and chaos as long as 6 > 1 . The difference between Eq. (|^) 
and the empirical quantum mechanical result is due to the approximation; 
also, such border lines are not expected to coincide exactly for the classical 
and quantum mechanical case. We may combine this relation with E ^ 
/(I — XI) {rir = 0,huj = 1, N and / sufficiently large) to obtain the estimate 
XI ~ 1/2 as a border line at around 6 = 1. Therefore, values of / and A which 
obey the condition 1/2 < A ■ / < 1 correspond to chaotic motion. (Note 
that, since outside the chaotic region the values of the angular momentum 
are almost sharp, values of / can be used to delineate different regions even 
though it is a dynamical variable.) As a consequence, the values used for A 
in nuclear physics, when applied for large metallic clusters, lead to a strong 
mixing of levels at the high end of the spectrum (see Fig. (7a, 7b)) even for 
values of b close to unity. Hence, shell structure is expected to be weakly 
pronounced for a large particle number even for mild deformation. For 6 = 2 
the emerging shell structure is confined to the very low end of the spectrum 
and to rather small values of A. 

The results found for the specific examples 6 = 1.1 and 6 = 2 do not 
change for negative parity states or for A > 0. The spectrum for larger 
values of A nicely slots in with the spectrum for A = to produce typical 
shell structure as long as b is close to unity. At the end of Section 3 we 
pointed out that a regular pattern begins to emerge for very large values of 
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p^p even when 6 = 2. Note also that, for large values of A (the quantum 
mechanical value of p^), (L)^ becomes diagonally dominant as seen from 
Eq.(|T2|) thus leading to lesser pronounced level repulsion in line with the 
increasingly regular classical situation. While levels of large A are submerged 
in the morass of a chaotic spectrum pertaining to lower A-values, they are 
expected to yield discernible shell structure in a Strutinsky type analysis. In 



calculations by [T^, ^, levels belonging to rather high shells occur below the 
Fermi surface and have to be taken into account in the Strutinsky averaging 
procedure; pure oscillator states (A = 0) are used for the upper part of the 
spectrum. Even for relatively small clusters with sizes up to A = 270 atoms 
[p3] the convergence criterium for the Strutinsky's shell correction method 



2^ required inclusion of spherical shell numbers up to A^ = 16. In essence, 
the occurrence of large A is associated with very high shell numbers, a region 
where the applicability of the model is no longer certain. For example, in 
[|T^ it was found that for large spherical clusters {A = 676 and 832) the 
energy minima are more pronounced than in the KS calculations, which is an 
artifact. In other words, while in the axial MHO potential the re-emergence 
of order and hence shell structure is a fact for large A, its physical impact 
cannot be assessed from this model alone. Moreover, even small non-axial 



deformations, which are discernible for small clusters [^, |2^, will destroy 
shell effects for large clusters, since L^ is no longer a constant of motion, i.e. 
A is no longer a good quantum number. 

5 Summary 

Using the classical analysis of single particle motion in the MHO potential our 
attention is focused on the distinction of regions in phase space of predomi- 
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nantly regular or chaotic motion, when the deformation b and the term (L)^ 
act simultaneously. Important similarities, not obvious at first glance, with 
regard to the shortest periodic orbits which define the gross shell structure of 
the single particle spectrum in the MHO potential and in the Woods-Saxon 
potential have been found previously ^l| . Typical trajectories of periodic 



orbits are displayed in Fig. (3). Both cases give rise to chaotic motion when 
the deformation in the Woods-Saxon and the (L)^-term in the MHO poten- 
tial are turned on. The chaotic behavior is of lesser significance, as long as the 
interest is focused only on the lower end of the spectrum and the deformation 
is fairly mild, that is for small metallic clusters. For large metallic clusters, 
however, with the much larger particle numbers, a higher range of the spec- 
trum becomes relevant; depending on the deformation, chaotic behavior may 
well interfere with the search for shell structure in this region. 

Regular and chaotic behavior are essentially separated by the curve EX ~ 
huj/{4b'^). When b > 1 and A = 0, the problem is regular and so is the vari- 
ation of the classical angular momentum; in fact the variation is determined 
by the frequencies u and u/b. If A is switched on, the variation becomes 
quickly irregular, but more importantly, the range of variation increases con- 
siderably with increasing A as it should be expected, since the term A(L)^ 
adds angular momentum. In quantum mechanics, switching on the (L)^ term 
not only lifts the degeneracies but brings about a strong mixing among all 
levels of the unperturbed basis. Therefore, the unavoidable truncation of the 
basis in such a situation must be handled with great circumspection to avoid 
the occurrence of spurious shell structure. While it is common wisdom that 
questions of truncation are delicate, we stress that the classical finding, which 
is the large fluctuation of the angular momenta, provide essential guidance 
to the problem of the truncation. 
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From Figs. 7 and from the relation EX ^ 1/(46^) {huj = 1) we conclude 
that shell structure is possible at higher energies (larger number of particles) 
only for accordingly small values of A. For A > 0.007 (which is smaller by one 
order of magnitude than the typical ones) and b = 2 hard chaos occurs already 
for E > 10 which is seen from the Wigner distribution of the nearest neighbor 
statistic. Shell structure is retrieved for lesser deformation at the same values 
of energy and A. One of the basic conclusions of our analysis is that using 
the MHO model for metallic clusters with a coupling strength A similar to 
that in nuclear physics would conflict, at least for larger deformations b, with 
observation of shell structure. Within the model considered we may speculate 
that large particle systems tend to restore the near-spherical symmetry to 
allow for shell structure which basically means increased stability for the 
appropriate particle number. There seem to be indications to this effect 
in calculations of the equilibrium shapes for sodium clusters with A < 300 
(see, for example |2^) using the MHO. We conclude that the MHO potential 



cannot yield shell structure for 6 > 2 and shell number larger than about 
4, or 6 ~ 1.1 and shell number larger than about 20. Energy minimization 
is therefore expected to yield decreasing deformation for increasing particle 
number. 

When the model was introduced originally in nuclear physics [|l^], nei- 



ther the unbounded motion nor the inherent chaotic nature was of interest. 
The (L)^ term has been introduced to modify the nuclear potential with an 
effective A-dependence. Its chaotic nature was masked at the time, since 
only A < 200 was of interest. The transfer of a well established nuclear 
structure model to the new physical situation of metallic clusters requires 
a re-assessment of the phenomenological model, since it entails systems of 
a considerably larger number of particles, where even a moderate deforma- 
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tion can have a pronounced effect. The connection between shell structure 
phenomena in the quantum mechanical spectrum and ordered motion in the 
analogous classical case leads to the conclusion that the model is useful but 
only in a restricted domain. Our analysis shows that the two aspects, the 
unboundedness and chaotic motion, are independent of each other. In fact, 
the regular classical motion for large values of (L)^ (/ > Zcrit = 1/-^) gives rise 
to quantum mechanical states which have virtually sharp values of / and are 
essentially orthogonal upon the chaotic states which, in turn, are character- 
ized by the inequality 1/2 < /A < 1. The other regular part of the spectrum 
which is displayed in Fig. 7 is obtained for IX < 1/2. 

The inclusion of higher multipoles into the potential for large system 
cannot change this finding. The chaotic nature inferred by adding the term 
(L)^ destroys the shell structure which is produced in a subtle way by the 
octupole deformation when added to the quadrupole deformed HO |3^. The 
hexadecapole deformation, when added to the quadrupole deformed HO, 
induces shell structure with a pattern similar to that of a lesser deformed 
quadrupole HO [^. And the combination yields shell structure for specific 



values of the octupole and hexadecapole strength [0. Since the (L)^ term 
dissolves all bunching of levels in a deformed potential at higher energies, 
i.e. for large particle numbers, shell structure can be restored only for a near 
spherical shape. This property is shared with the effect of the hexadecapole 
term mentioned above. In fact, similarities between the two terms have 
been pointed out in [|13|. Other shapes of potentials [^ seem to indicate 



that adding higher multipoles gives rise to shell structure for large particle 
numbers. It is a challenge to understand such claims in terms of periodic 
orbits. Work towards this aim is in progress. 
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n 


En 


L2 


AL2 


1 


1.25 


0.25 


1.35 


2 


2.25 


5.25 


5.73 


3 


3.25 


10.3 


11.5 


4 


3.25 


2.75 


4.85 


5 


4.25 


17.8 


15.7 


6 


4.25 


15.3 


17.3 


7 


5.25 


32.8 


29.1 


8 


5.25 


20.3 


23.1 


9 


5.25 


5.35 


8.78 


17 


8.25 


80.3 


70.7 


18 


8.25 


77.8 


69.9 


19 


8.25 


42.7 


38.2 


20 


8.25 


35.3 


40.6 


26 


10.25 


130 


115 


27 


10.25 


113 


99 


28 


10.25 


108 


97 


29 


10.25 


56 


50 


30 


10.25 


45 


52 



Table (1). Average values and variances of the operator {Vf for a few low 
energy states for 6 = 2 and A = 0. 
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n 


En 


L2 


AL^ 


1 


1.25 


0.29 


1.47 


2 


2.19 


5.85 


6.51 


3 


3.12 


13.3 


14.7 


4 


3.22 


2.82 


5.24 


5 


4.04 


24.0 


28.9 


6 


4.06 


18.7 


16.3 


7 


4.85 


43.8 


47.3 


8 


4.96 


37.0 


39.7 


9 


5.19 


6.41 


9.80 


13 


5.91 


32.3 


36.6 


14 


6.25 


1657 


1862 


15 


6.30 


2163 


1973 


16 


6.50 


1947 


1963 


27 


7.76 


51.3 


117 


28 


7.93 


1504 


1740 


29 


7.95 


1413 


1810 


30 


8.13 


836 


1479 



Table (2). Same as Table (1) but for A = 0.01. 
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n 


En 


L2 


(0 


AL2 


1 


-28.1 


506 


(22) 


9.6 


2 


-21.5 


420 


(20) 


9.0 


3 


-19.5 


420 


(20) 


9.9 


4 


-15.6 


342 


(18) 


8.4 


9 


-6.7 


272 


(16) 


9.7 


10 


-6.2 


210 


(14) 


7.4 


17 


-0.17 


110 


(10) 


6.7 


18 


1.1 


156 


(12) 


8.7 


20 


1.45 


0.01 


(0) 


0.2 


21 


1.87 


71.6 


(8) 


6.9 


22 


1.89 


110 


(10) 


7.6 


23 


2.7 


6.7 


(2) 


4.4 


24 


2.9 


39 


(6) 


10.5 


25 


3.0 


156 


(12) 


9.7 


26 


3.3 


22.8 


(4) 


8.5 


27 


3.4 


0.2 


(0) 


1.7 



Table (3). Same as Table (1) but for b = 1.1 and A = 0.1. Note that for 
this small deformation the levels are nearly good eigenstates of (L)^ with 
averages close to /(/ + 1) and small variances. 
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Figure Captions 

Fig.l Typical simple periodic orbits for b = 1. The orbits with the loops 
at the corners are for A > huj/{4E) while the other two are for A = huj / {25E) . 

Fig. 2 Surfaces of section in the plane g = 0. In (a) the dependence on 
A of accessible phase space (shaded) is illustrated; left: A < Tiu /{AE), right: 
A > huj/{AE). In (b) four orbits of different initial conditions are displayed 
for A < hiu/{AE). 

Fig. 3 Typical short periodic orbits in the g — z— plane. Reference to Fig. 2 
is made in the main text. Figs.(c) and (d) are self-tracing orbits. 

Fig. 4 Two orbits for the same choice of parameters (6 = 2, A = 0.01). 
The left hand column displays the surface of section and the corresponding 
orbit for a regular orbit while the right hand column gives the respective 
plots for a typical chaotic orbit. 

Fig. 5 Two orbits for the same choice of parameters and initial conditions 
but for different values of p^. The left hand column displays the surface of 
section at g = 8 and the corresponding orbit for a regular orbit at p^p = 61 
while the right hand column gives the respective plots for a typical chaotic 
orbit at j9(p = 1. 

Fig. 6 A nearest neighbour distribution of the unfolded spectrum for 6 = 2 
and A = 0.01. 

Fig. 7 (a) Lower end (about 6% of all levels) of the energy spectrum in 
units of huj as a function of the strength A for states whose eigenvectors 
exhaust 50% of the norm in the first fifth of the total length. The dividing 
line EX = 1/(46^) can be nicely discerned, (b) Lower end of the spectrum 
for all levels, (c) A section of (a), (d) Same as (c) but only for states whose 
eigenvectors exhaust 80% of the norm in the first fifth of the total length. 



30 



(a) 



P 



^ri 


hP 


J^ 


WV 




(b) 

P 

10k ^. 



-6 ■. l-A -2 



^ Z 



-lO' 



-10 



20 




^ Ul 




10 




5 




-5 




-10 




-lb" 




-20 





\ 



10 



. 


p 

20r 

15- 


•! / 


•J 


.•"•10. 

•••i 


•k 
_ •.^ 




-20A -10 


•■ -5 


• .* 
« 


10 I 20 
^^'^^ •' 


■ ■ •• 


-20 




• ." . 







7.5 

7 

6.5 


1 .'^ 

/ / 

/ 
/ 


\ 


-1 

5.5 


v 




.J r. :■:■':. .20 
: ■■.•.>-'" ••-■••• .• 






,v'-'*.^-. fa ..••■•• 1.5 V 


• * .:• ."2 







0.2 0.4 0.6 0, 
100 • strength 





(b) 




~-~=^^^^^ ■.■'■■ ^rr^»^:^£;itil?Ct5:; 


14 


— ==^ss8 ■ .'■ ■.-•>:.--vv,.'^T^/;*-'Vvv?^<:^i5 






12 








10 


■ ■ — ■ , ■ v,,.\.' 




— ^^■'Wt^^fw^^® 


5 


^W\tttWm^ 




~\ V\\Vi\"\t^ttl?^^^^^^ 


6 


) AAvT\V\VwX^;W#;^^^^m^ 




V V \W^\W\?^W^^S^^^^ 


4 


\ \ \\\\\\\\\\\\\\\\C\\\\\^ 




\ \ \\\\\\\\\\\%W\W 


2 


\ \ v\W{\^M^s!t^^^ " 



0.2 0.4 0.6 
100 • strength 



15 
14 




(c) 








"■''■i!;;;:;;:;;;!;!' -■ 


■"''?-:'>:'5::;;:;:7!K:w!;!J^^^^^^ 




;::5,:'!;,-;i:;::'|j;i.i 


3i . 
III . 


>1 1 o 

U 
Q) 

% 12 


'*«iii;;;:£=:::::::::::; 


^SsiSSi: 


i:~;:::::|::ft 
Jill;: 




11 
1 n 


^^^^Sigg::::::: 


*s:::*:«:;;:;;^''''''-%?;'i; 




t'-'iSiilSJ; 


i*::;:*;;i!:: 



0.2 0.4 0.6 0.8 1 
100 • strength 



15 



14- 



>1 1 o 

u 

Q) 

S 12 



11 



10 



iig-aa::::'" 



(d) 









0.2 0.4 0.6 0.8 1 
100 • strength 



